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Abstract 

'Constituent quarks' means massive quarks, in contrast to the nearly massless u, d, s quarks of 
the QCD lagrangian. The dynamical or constituent masses appear owing to the spontaneous chiral 
symmetry breaking in QCD. A realistic mechanism of the chiral symmetry breaking is provided by 
instantons. Therefore, I review the present status of the QCD instanton vacuum, emphasizing the 
^ ■ mechanism of chiral symmetry breaking and of the generation of the constituent quark mass. I end 
up with the Nambu-Jona-Lasinio-type effective theory to which QCD is reduced at low momenta. 
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1 Introduction 



Quantum Chromodynamics (QCD) is the theory of strong interactions. From a theorist's point of 
view strong interactions means that all dimensionless quantitites of the theory are, generally speaking, 
of the order of unity. Meanwhile, there are at least two fundamental facts about strong interactions 
showing that actually it is not exactly so: there are still certain small parameters hidden. One is 
that nuclei can be, to a good accuracy, viewed as made of weakly bound nucleons — in a true 'strong 
interactions' case they would rather look like a quark soup. The second puzzle is that the nucleon 
itself can be, to a good accuracy, viewed as built of three constituent quarks — in a true 'strong 
interactions' case it would rather look like a pack of an indefinite number of quarks, antiquarks and 
gluons. 

The second puzzle (and probably the first, too) is qualitatively explained by the spontaneous breaking 
of the (approximate) chiral symmetry of QCD. The bare ot current masses of the light quarks 
{m u ~ 4 MeV, ~ 7 MeV, m s ~ 150 MeV) are small as compared to the typical scale of strong 



interactions, say, the p-meson mass, 770 MeV, and can be put to zero. In that limit (called the 
chiral limit) the QCD lagrangian possesses an 18-parameter symmetry, called chiral symmetry It is 
the symmetry under independent U(3) x U(3) rotations of the left- and right-handed components of 
the u, d, s fields. Equivalently, it is the symmetry under U{3)y ector x U (3) Axial rotations. Were this 
symmetry exact, we would observe degeneracy between states with opposite parity but otherwise 
the same quantum numbers. For example, the vector p-meson would be degenerate with the axial 
ai-meson; the nucleon (| + , 940) would be degenerate with the resonance (| , 1535), etc. Since it is 
not the case, we conclude that the U{3)a symmetry is broken spontaneously. The order parameter 
for symmetry breaking is the quark or chiral condensate, 

tyf) ~ -(250 MeV) 3 , ip = u,d, s. (1) 
Its U(3)a rotations produce nine massless Goldstone bosons - the pseudoscalar nonet. 

It should be noted that the above quantity is well defined only for massless quarks, otherwise it 
is somewhat ambigious. By definition, this is the quark Green function taken at one point; in 
momentum space it is a closed quark loop. Had the quark propagator only the "slash" term, the 
trace over the spinor indices understood in this loop would give an identical zero. Therefore, chiral 
symmetry breaking implies that a massless (or nearly massless) quark develops a non-zero dynamical 
mass (i.e. a "non-slash" term in the propagator). There are no reasons for this quantity to be a 
constant independent of the momentum; moreover, we understand that the dynamical mass should 
anyhow vanish at large momenta. Its value at zero momentum can be estimated as one half of the 
p meson mass or one third of the nucleon mass, that is about 

M(0) ~ 350 - 400 MeV. (2) 

This scale is also related to chiral symmetry breaking and should emerge together with the conden- 
sate. 

The spontaneous breaking of chiral symmetry is the main dynamical happening in QCD: it determines 
the face of the strong interactions world. Indeed, it explains why the pseudoscalar mesons are light 
(they are called pseudo-Goldstone bosons since they get a non-zero mass from the small current 
quark masses which break somewhat the chiral symmetry from the beginning), and why nucleons 
are heavy and not degenerate in parity. The appearance of the dynamical quark mass (I prefer this 
word to the constituent mass), related to the chiral symmetry breaking, justifies, qualitatively, the 
old non-relativistic or semi-relativistic models of hadrons. Indeed, once one gets @, one can state 
that mesons are predominantly made of a quark and an antiquark, and baryons are made of three 
quarks: an addition of a quark- antiquark pair costs 700 — 800 MeV which is a lot. Moreover, the 
lightest vector, axial and tensor mesons and baryons should be rather loosely bound: first, because 
their masses are not far away from the sum of the dynamical quark masses, second, because, after 
the appearance of a rather large scale in the theory, there is no place for the "infrared catastrophy" 
- the gluon coupling constant a s remains finite and in fact, as known from hard hadronic processes, 



relatively small. Therefore, one would expect that neither the one-gluon exchange nor the linear 
confining potential plays a crucial role in determining the world of light hadrons (I shall soon cite 
direct evidence from lattice experiments substantiating this view). Of course, one would need explicit 
confining forces to describe highly excited resonances lying on linear Regge trajectories, when hadrons 
are artificially "stretched" by centrifugal forces, but they are not so much important for describing 
low-lying hadrons. The dominant phenomenon is chiral symmetry breaking, and if we want to 
understand strong interactions we have to understand how to get ([l], |2|) from the only dimensional 
quantity in QCD, that is Aqcd which is the ultraviolet cut-off-independent combination of the cut-off 
fi and the gauge coupling constant g 2 (n) given at that cut-off: 

11 2 

b = Y Nc ~ z Nf > (3) 

where N c is the number of colours (N c = 3) and Nf is the number of active quark flavours. The value 
of Aqcd depends on the regularization scheme used; in the MS scheme Ajjg ~ 250 — 300 MeV. 

There is a remarkable puzzle, called the U(l) paradox (Weinberg, 1974): eight out of nine pseu- 
doscalar mesons (tc, K, rf) look indeed like pseudo-Goldstone bosons and obey the PCAC relations 
while the singlet ?/-meson does not (the PCAC predicts its mass to be about 400 MeV vs. 958 MeV 
(exper.)). A related puzzle is the ?7 — ^ 37r decay. It has been known from 't Hooft (1976) that, qual- 
itatively, one needs instantons to cure the U(l) paradox (see also (Diakonov, 1982)). Probably the 
modest 7/-meson tells us about strong interactions more than any other hadron! It says that one 
has to take into account instanton configurations in QCD: there is no other way to get a heavy rj'- 
meson. We shall see that instantons bring in a most realistic scenario of spontaneous chiral symmetry 
breaking as well. 

2 Instantons 

Instantons are certain configurations of the Yang-Mills potentials for the gluon field A a (x), satisfying 
the equations of motion D^F^ = in euclidean space, i.e. in imaginary time. The solution has been 
found by Belavin, Polyakov, Schwartz and Tiupkin (1975); the name "instanton" has been suggested 
by 't Hooft in (1976), who also made a major contribution to the investigation of the instantons 
properties. 

In QCD instantons are the best studied non-perturbative effects, leading to the formation of the gluon 
condensate (Shifman, Vainshtein and Zakharov, 1979) and of the so-called topological susceptibility 
needed to cure the U(l) paradox. The QCD instanton vacuum has been studied starting from the 
pioneering works in the end of the seventies (Callan, Dashen and Gross, 1978; Shuryak, 1982); a 
quantitative treatment of the instanton ensemble, based on the Feynman variational principle, has 
been developed by Diakonov and Petrov (1984a). The most striking success of the QCD instanton 
vacuum is its capacity to provide a beautiful mechanism of the spontaneous chiral symmetry breaking 
(Diakonov and Petrov, 1984b, 1986a). Moreover, the instanton vacuum leads to a very reasonable 
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effective chiral lagrangian at low energies, including the Wess-Zumino term, etc., which, in its turn 
gives a nice description of nucleons as chiral quark solitons (Diakonov and Petrov, 1986b). 

A detailed numerical study of dozens of correlation functions in the instanton medium undertaken 
by Shuryak and collaborators (1993a) (earlier certain correlation functions were computed analyt- 
ically by Diakonov and Petrov (1984b, 1986b)) demonstrated an impressing agreement with the 
phenomenology (Shuryak, 1993b) and, more recently, with direct lattice measurements (Chu et al, 
1994). As to baryons, the instanton-motivated chiral quark soliton model (Diakonov and Petrov, 
1986b) also leads to a very reasonable description of dozens of baryon characteristics (see Klaus 
Goeke's contribution to these Proceedings). 

More recently the instanton vacuum was studied in direct lattice experiments by the so-called cooling 
procedure (Teper, 1985; Ilgenfritz et al, 1986; Polikarpov and Veselov, 1988; Chu et al, 1994). It was 
demonstrated that instantons and antiinstantons (Ps and Ps for short) are the only non-perturbative 
gluon configurations surviving after a sufficient smearing of the quantum gluon fluctuations. The 
measured properties of the II medium appeared (Chu et al, 1994; Michael and Spencer, 1995) to 
be close to that computed from the variational principle (Diakonov and Petrov, 1984a) and to what 
had been suggested by Shuryak (1982) from phenomeno logical considerations. 

Cooling down the quantum fluctuations above instantons kills both the one-gluon exchange and the 
linear confining potential (a small residual string tension observed in the cooled vacuum (Polikarpov 
and Veselov, 1988; Chu et al, 1994) is probably due to the instanton-induced rising potential at 
intermediate distances (Diakonov, Petrov and Pobylitsa, 1989). Nevertheless, in the cooled vacuum 
where only Ps and Ps are left, the correlation functions of various mesonic and baryonic currents, as 
well as the density-density correlation functions representing the quark wave functions in hadrons, 
appear to be quite similar to those of the true or "hot" vacuum (Chu et al, 1994). 

In recent lectures (Diakonov, 1995) I have reviewed the building blocks of the instanton vacuum; an 
interested reader is addressed to those lectures or original literature. Here I shall review the main 
facts about instantons, in particular the way one gets chiral symmetry breaking and the 'constituent' 
quarks from instantons. 



3 Instanton field 

Physically, one can think of instantons in two ways: on one hand it is a tunneling process occuring in 
imaginary or euclidean time (this interpretation belongs to V.Gribov, 1976), on the other hand it is 
a localized pseudoparticle in the euclidean space (A.Polyakov, 1977). The gluon field of the instanton 
in the singular gauge is ('t Hooft, 1976): 

2o 2 O ab fi i (r — z) 

A{ a = , P N 'r V , ^— , O ab = Tr(UH a Ua l ), O ai O aj = 5 ij . (4) 

M (x — z) 2 [p 2 + (x — z) l \ 



Here the 4-dim. vector z^ is called the instanton centre, p is the instanton size, the rectangular 
matrix O a \ a — 1, , , (N% — 1), i — 1, 2, 3) gives the orientation of the instanton in the colour SU(N C ) 
space. All in all there are 



4 (centre) + 1 (size) + (4iV c — 5) (orientations) = 4N C (5) 

so called collective coordinates desribing the field of the instanton, of which the action is independent. 
The tensors ff^ are called 't Hooft symbols ('t Hooft, 1976). 

The field strength of an instanton (centered at z^ = 0) is 

F a - 4p2 Q ai (rf - Irf XaXv - Irf X ^ Xf3 \ F a F a - 192p4 (R\ 

and satisfies the self-duality equation F" = F° . The anti-instanton satisfies the antz-self-dual 
equation, F = —F; it is given by eqs.(|4], ^|) with the replacement fj — > rj. 

The action of one (anti) instanton is 



S=- I d 4 x F^F% = Stt 2 . (7) 



4 Gluon condensate 



The QCD perturbation theory implies that the fields A^(x) are performing quantum zero-point 
oscillations; in the lowest order these are just plane waves with arbitrary frequences. The aggregate 
energy of these zero-point oscillations, (B 2 + E 2 )/2, is divergent as the fourth power of the cutoff 
frequency, however for any state one has (F^ u ) = 2(B 2 — E 2 ) = 0, which is just a manifestation 
of the virial theorem for harmonic oscillators: the average potential energy is equal to that of the 
kinetic (I am temporarily in the Minkowski space). One can prove that this is also true in any 
order of the perturbation theory in the coupling constant, provided one does not violate the Lorentz 
symmetry and the renormalization properties of the theory. Meanwhile, we know from the QCD 
sum rules phenomenology that the QCD vacuum posseses what is called gluon condensate (Shifman, 
Vainshtein and Zakharov, 1979): 



32^(^ a ,> - (200 MeVf > 0. (8) 

Instantons suggest an immediate explanation of this basic property of QCD. Indeed, instanton is a 
tunneling process, it occurs in imaginary time; therefore in Minkowski space one has Ef = ±iB? 
(this is actually the self-duality equation). In euclidean space the electric field is real as well as the 
magnetic one, and the gluon condensate is just the average action density. Let us make a quick 
estimate of its value. 



Let the total number of 7's and J's in the 4-dimensional volume V be N. Assuming that the average 
separations of instantons are larger than their average sizes (to be justified below), we can estimate 
the total action of the ensemble as the sum of invidual actions: 



WV = J d 4 xF* u ^N-32ir 2 , (9) 

hence the gluon condensate is directly related to the instanton density in the 4-dimensional euclidean 
space-time: 

1 N 1 

(F a F a ) = — (10) 

32tt 2 V 1X1/1 ~ V ~ R 4 1 ; 

In order to get the phenomenological value of the condensate one needs thus to have the average 
separation between pseudoparticles (Shifman, Vainshtein and Zakharov, 1979; Shuryak, 1982) 



R ~ = 1 fm. (11) 

200 MeV J V ; 

There is another point of view on the gluon condensate which I describe briefly. In principle, all 
information about field theory is contained in the partition function being the functional integral 
over the fields. In the euclidean formulation it is 



Z = j DA.exp (-± J d*xF^ e~ ET , (12) 

where I have used that at large (euclidean) time T the partition function picks up the ground state 
or vacuum energy E. For the sake of brevity I do not write the gauge fixing and Faddeev-Popov 
ghost terms. If the state is homogeneous, the energy can be written as E = OuV^ where 9^ is the 
stress-energy tensor and is the 3-volume of the system. Hence, at large 4-volumes V = V^T 
the partition function is Z = exp(— 644V). This 644 includes zero-point oscillations and diverges 
badly. A more reasonable quantity is the partition function, normalized to the partition function 
understood as a perturbative expansion about the zero-field vacuum. The latter can be distinguished 
from the former by imposing a condition that it does not contain integration over singular Yang-Mills 
potentials: instanton potentials are singular at the origins. One has 



(13) 

We expect that the non-perturbative vacuum energy density 6*44 — 9^ is a negative quantity since 
we have allowed for tunneling: as usual in quantum mechanics, it lowers the ground state energy. If 
the vacuum is isotropical, one has 6*44 = 0^/4. Using the trace anomaly, 
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~ —b 

one gets (Diakonov and Petrov, 1984a): 



6 W ~ -b-^r, b = —N c , (14) 



= exp \^V(F 2 j32n 2 ) NP j (15) 

where {F 2 v )np is the gluon field vacuum expectation value which is due to non-perturbative fluctu- 
ations, i.e. the gluon condensate. The aim of any QCD-vacuum builder is to minimize the vacuum 
energy or, equivalently, to maximize the gluon condensate. 



5 One-instanton weight 



The words 'instanton vacuum' mean that one assumes that the QCD partition function is mainly 
saturated by an ensemble of interacting J's and J's together with quantum fluctuations about them. 
Instantons are necessarily present in the QCD vacuum if only because they lower the vacuum energy 
in respect to the purely perturbative (divergent) one. The question is whether they give the dominant 
contribution to the gluon condensate, and to other basic quantities. To answer this question one has 
to compute the partition function ([12]) assuming that it is mainly saturated by instantons, and to 
compare the obtained gluon condensate with the phenomenological one. This work has been done a 
decade ago in ref. (Diakonov and Petrov, 1984a); today direct lattice measurements seem to confirm 
that the answer to the question is positive: the observed density of J's and J's is in agreement with 
the estimate ([TT|). 

The starting point of this calculation is the contribution of one isolated instanton to the partition 
function (0), or the one-instanton weight. To get a reasonable result it must be i) normalized 
(to the determinant of the free quadratic form, i.e. with no background field), ii) regularized (for 
example by using the Pauli-Villars method), and ra) accounted for the zero modes. The resulting 
one-instanton contribution to the partition function (normalized to the free one) is an integral over 
the 4N C collective coordinates of an instanton ('t Hooft, 1976; Bernard, 1979): 
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(16) 
(17) 



The product of the last two factors in eq. (|T6| ) is actually a combination of the cut-off p and the 
bare coupling constant g 2 (p) given at this cut-off, which is cutoff-independent; it can be replaced by 
(pAqcd) UNc ^ 3 , see eq. (|j). This is the way Aqcd enters into the game; henceforth all dimensional 
quantities will be expressed through Aqcd, which is, of course, a welcome message. The numerical 
coefficient C(N C ) depends explicitly on the number of colours; it also implicitly depends on the 
regularization scheme used. 

Note that the g 2 in the pre-exponent starts to "run" only at the 2-loop level, hence its argument 
is taken at the ultra-violet cut-off p. The 2-loop instanton weight can be found in refs. (Diakonov, 
1995; Diakonov, Polyakov and Weiss, 1955). 



In both one- and two-loop approximations the integral over the instanton sizes p in eq. ( |T7D diverges 
as a high power of p at large p: this is of course the consequence of asymptotic freedom. It means 
that individual instantons tend to swell. This circumstance plagued the instanton calculus for many 
years. If one attemts to cut the p integrals "by hand", one violates the renormalization properties of 
the YM theory. Actually the size integrals appear to be cut from above due to instanton interactions. 



6 Instanton ensemble 



To get a volume effect from instantons one needs to consider an // ensemble, with their total number 
./V proportional to the 4-dimensional volume V. Immediately a mathematical difficulty arises: any 
superposition of J's and J's is not, strictly speaking, a solution of the equation of motion, therefore, 
one cannot directly use the semiclassical approach of the previous section. There are two ways to 
overcome this difficulty. One is to use a variational principle (Diakonov and Petrov, 1984a), the other 
is to use the effective YM lagrangian in the instanton field (Diakonov and Polyakov, 1993). 

The idea of the variational principle is to use a modified YM action for which a chosen II ansatz is 
a saddle point. Exploiting the convexity of the exponent one can prove that the true vacuum energy 
is less than that obtained from the modified action. One can therefore use variational parameters 
(or even functions) to get a best upper bound for the vacuum energy. We call it the Feynman 
variational principle since the method was suggested by Feynman in his famous study of the polaron 
problem. Todays direct lattice investigation of the II ensemble seem to indicate that Petrov and I 
have obtained rather accurate numbers in this terrible problem. Therefore I will cite the numerics 
from those calculations in what follows. 

The main finding is that the II ensemble stabilizes at a certain density related to the Aqcd parameter 
(there is no other dimensional quantity in the theory!): 



(^/32tt 2 ) ~ - > (0.75A^) 4 (18) 

which would require ~ 270 MeV to get the phenomenological value of the condensate. It 

should be mentioned however that using more sophisticated variational Ansatze one can obtain a 
larger coefficient in eq. (|I8|) and hence would need smaller values of A. 

The average sizes p appears to be much less than the average separation R. Numerically we have 
found for the SU(3) colour: 



which coincides with what was suggested previously by Shuryak (1982) from considering the phe- 
nomenological applications of the instanton vacuum. This value should be compared with that found 
from direct lattice measurements (Chu et al, 1994): p/R ~ .37 — .4, depending on where one stops 
the cooling procedure. The packing fraction, i.e. the fraction of the 4-dimensional volume occupied 



by instantons apears thus to be rather small, 7r 2 p 4 /-R 4 ~ 1/8. This small number can be traced back 
to the "accidentally" large numbers appearing in the 4-dimensional YM theory: the lliV c /3 of the 
Gell-Mann-Low beta function and the number of zero modes being 4N C . Meanwhile, it is exactly 
this small packing fraction of the instanton vacuum which gives an a posteriori justification for the 
use of the semi-classical methods. As I shall show in the next sections, it also enables one to identify 
adequate degrees of freedom to describe the low-energy QCD. 



7 Chiral symmetry breaking by instantons: qualitative deri- 
vation 

The key observation is that the Dirac operator in the background field of one (anti) instanton has an 
exact zero mode with A = ('t Hooft, 1976). It is a consequence of the general Atiah-Singer index 
theorem; in our case it is guaranteed by the unit Pontryagin index or the topological charge of the 
instanton field. These zero modes are 2-component Weyl spinors: right-handed for instantons and 
left-handed for antiinstantons. 

For infinitely separated I and 7 one has thus two degenerate states with exactly zero eigenvalues. As 
usual in quantum mechanics, this degeneracy is lifted through the diagonalization of the hamiltonian, 
in this case the hamiltonian is the full Dirac operator. The two "wave functions" which diagonalize 
the "hamiltonian" are the sum and the difference of the would-be zero modes, one of which is a 
2-component left-handed spinor $i, and the other is a 2-component right-handed spinor <3>2- The 
resulting wave functions are 4-component Dirac spinors; one can be obtained from another by mul- 
tiplying by the 75 matrix. As the result the two would-be zero eigenstates are split symmetrically 
into two 4-component Dirac states with non-zero eigenvalues equal to the overlap integral between 
the original states $i ;2 : 

A = ±|T 12 |, T 12 = / fl)Q 2 R ^°° - 2 -^PlTr{U\U 2 Rt 2 ). (20) 

We see that the splitting between the would-be zero modes fall off as the third power of the distance 
between I and 7; it also depends on their relative orientation. 

When one adds more Ps and 7's each of them brings in a would-be zero mode. After the diagonal- 
ization they get split symmetrically in respect to the A = axis. Eventually, for an II ensemble 
one gets a continuous band spectrum with a spectral density u(X) which is even in A and finite at 
A = 0. Meanwhile, there is a general formula relating the chiral condensate with the average spectral 
density of the Dirac operator at zero A (Banks and Casher, 1980; Diakonov and Petrov, 1984b): 

(W) = -^sign(m) tt^O) (21) 
Since Ps and 7's lead to a non-zero z/(0) they automatically give the non-zero chiral condensate, 



which signals chiral symmetry breaking. 

One can make a quick estimate of (ipip): Let the total number of I's and J's in the 4- dimensional 
volume V be N. The spread A of the band spectrum of the would-be zero modes is given by their 
average overlap (p0|): 



A ~ (d*R/V)T(R)T*(R) ~ (22) 

where p is the average size and R = (N/V)~ l l 4 is the average separation of instantons. Note that 
the spread of the would-be zero modes is parametrically much less than 1/p which is the typical 
scale for the non-zero modes. Therefore, neglecting the influence of the non-zero modes is justified if 
the packing fraction of instantons is small enough. From eq. (pl"D one gets an estimate for the chiral 
condensate induced by instantons: 



(U) = ->) - -f £ ~ -i. (23) 

It is amusing that the physics of the spontaneous breaking of chiral symmetry resembles the so- 
called Anderson conductivity in disordered systems. Imagine random impurities (atoms) spread over 
a sample with final density, such that each atom has a localized bound state for an electron. Due to 
the overlap of those localized electron states belonging to individual atoms, the levels are split into a 
band, and the electrons become delocalized. That means conductivity of the sample. In our case the 
localized zero quark modes of individual instantons randomly spread over the volume get delocalized 
due to their overlap, which means chiral symmetry breaking. 

I should mention that the idea that instantons can break chiral symmetry has been discussed previ- 
ously (see refs. (Caldi, 1977; Carlitz and Creamer, 1978; Callan, Dashen and Gross, 1978, Shuryak, 
1 982) , however the present mechanism and a consistent formalism has been suggested and developed 
in the papers (Diakonov and Petrov, 1984b, 1986a). 



8 Chiral symmetry breaking by instantons: quark propa- 
gator 

Having explained the physical mechanism of chiral symmetry breaking as due to the derealization of 
the would-be zero fermion modes in the field of individual instantons, I shall indicate how to compute 
observables in the instanton vacuum. The main quantity is the quark propagator in the instanton 
vacuum, averaged over the instanton ensemble. This quantity has been calculated in (Diakonov 
and Petrov, 1984b; Pobylitsa, 1989). In particular, Pobylitsa has derived a closed equation for 
the averaged quark propagator, which can be solved as a series expansion in the formal parameter 
Np 4 /VN C which numerically is something like 1/250. 



The result of these works is that in the leading order in the above parameter the quark propagator 
has the form of a massive propagator with a momentum- dependent dynamical mass: 



S(p) = ^rw)' M ^) = c Vi7ivr FM ' (24) 

where -F(^) is a combination of the modified Bessel functions and is related to the Fourier transform 
of the quark zero mode: it is equal to one at z = and decreases rapidly with the momentum, 
measured in units of the inverse average size of instantons; c is a constant of the order of unity which 
depends slightly on the approximation used in deriving the propagator. Note that the dynamical 
quark mass is non-analytical in the instanton density. 

Fixing the average density by the empirical gluon condensate (see section 4) so that R ~ 1 fm and 
fixing the ratio p/R — 1/3 from our variational estimate, we get the value of the dynamical mass at 
zero momentum, 



M (0) ~ (350) MeV (25) 



while the quark condensate is 



<*> = 1 / $F*v> = ~ iN < I $-<7Tmi) H -« 255 MeV)Z - (26) 

Both numbers, ( p5[) and (|26|), appear to be close to their phenomenological values. 



Using the above small parameter one can also compute more complicated quantities like 2- or 3-point 
mesonic correlation functions of the type 



(J A (x)J B (y)), (J A (x)J B (y)J c (z)), J A = i>T A il> (27) 

where T A is a unit matrix in colour but an arbitrary matrix in flavour and spin. Instantons influence 
the correlation function in two ways: i) the quark and antiquark propagators get dressed and obtain 
the dynamical mass, as in eq. fl24]), ii) quark and antiquark may scatter simultaneously on the same 
pseudoparticle; that leads to certain effective quark interactions. These interactions are strongly 
dependent on the quark-antiquark quantum numbers: they are strong and attractive in the scalar 
and especially in the pseudoscalar and the axial channels, and rather weak in the vector and tensor 
channels. 

Since we have already obtained chiral symmetry breaking by studying a single quark propagator in 
the instanton vacuum, we are doomed to have a massless Goldstone pion in the psuedoscalar and 
axial correlators. Having a concrete dynamical realization of chiral symmetry breaking at hand, we 
can not only check the general Ward identities of the PCAC (which work of course) but we are in 
a position to find quantities whose values do not follow from general relations. One of the most 
important quantities is the F n constant: it can be calculated as the residue of the pion pole. We get: 



const ( p x 
In = — r 
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In -~ 100 MeV vs. 93 MeV (exper.). (28) 
P 



JZJ \ 

This is a very instructive formula. The point is, F n is surprisingly small in the strong interactions 
scale which, in the instanton vacuum, is given by the average size of pseudoparticles, l/p — 600 MeV. 
The above formula says that F w is down by the packing fraction factor (p/R) 2 ~ 1/9. It can be said 
that F n measures the diluteness of the instanton vacuum! However it would be wrong to say that 
instantons are in a dilute gas phase - the interactions are crucial to stabilize the medium and to 
support the known renormalization properties of the theory, therefore they are rather in a liquid 
phase, however dilute it may turn to be. 

By calculating three-point correlation functions in the instanton vacuum we are able to determine, 
e.g. the charge radius of the Goldstone excitation: 



' 7T 



s/N c 



2?rF, 



(340 MeV)- 1 vs. (310 MeV)- 1 (exper.). (29) 



n 



A systematic numerical study of various correlation functions in the instanton vacuum has been 
performed by Shuryak and collaborators (1993a). In all cases considered the results agree well or 
very well with experiments and phenomenology. As I already mentioned in the introduction, similar 
conclusions have been recently obtained from direct lattice measurements (Chu et al., 1994). 



9 Chiral symmetry breaking by instantons: Nambu— Jona- 
Lasinio model 

The idea of the first two derivations of chiral symmetry breaking by instantons, presented above, is: 
"Calculate quark observables in a given background gluon field, then average over the ensemble of 
fields", in our case the ensemble of J's and 7's . The idea of the third derivation is the opposite: 
"First average over the II ensemble and obtain an effective theory written in terms of interacting 
quarks only. Then compute observables from this effective theory". This approach is in a sense 
more economical; it has been developed in (Diakonov and Petrov, 1986a; Diakonov, 1995; Diakonov, 
Polyakov and Weiss, 1995). 

Quark interaction arises when two or more quarks happen to scatter on the same pseudoparticle; 
averaging over its positions and orientations results in a four- (or more) fermion interaction term 
whose range is that of the average size of instantons. The most essential way how instantons influence 
quarks is, of course, via the zero modes. Since each massless quark flavour has its own zero mode, 
it means that the effective quark interactions will be actually 2Nf fermion ones. They are usually 
referred to as 't Hooft interactions as he was the first to point out the quantum numbers of these 
effective instanton-induced interactions ('t Hooft, 1976). In case of two flavours they are four- 
fermion interactions, and the resulting low-energy theory resembles the old Vaks-Larkin-Nambu- 



Jona-Lasinio model (Vaks et ah, 1961) which is known to lead to chiral symmetry breaking. 

Let me quote the result for the case of two flavours, Nf = 2 (Diakonov and Petrov, 1986a). In that 
case one gets a 4-quark interaction vertex: 

2N 2 r d 4 k 1 d 4 k 2 d%d% 
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where is the left-handed component of the quark of flavour / (/ = u,d). For the antiinstanton- 
induced vertex Y^~' one has to replace left handed components by right-handed. M(k) is the same 
dynamical quark mass as obtained in another approach outlined in the previous section. It is given 
by {z = hp/2): 

M(k) = M(0)F 2 (k), M(0) = A(2vrp) 2 , 
F(k) = 2z[I (z)K 1 (z) - h(z)K (z) - -J^K^z) k ^ F(0) = 1. (31) 

The value M(0) (or A) is found from the equation (Diakonov and Petrov, 1984b, 1986a) (called 
sometimes self-consistency or gap equation): 

4N C r d A k M 2 (k) . . 

c 1 1. (32) 
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It is seen from eq. ([32] ) that the dynamically generated mass is of the order of M(0) ~ \]N/ (VN c )p. 
Knowing the form of M(k) given by eq. (|3~1~D and using the "standard" values N/V = (1 fm) -4 , p = 
(1/3) fm we find numerically M(0) ~ 350 MeV. If one neglects M 2 in the denominator of eq. ( |3"2"D 
(Pobylitsa, 1989) one gets M(0) ~ 420 MeV. This deviation indicates the accuracy of the "zero mode 
approximation" used in this derivation: it is about 15%. 

Note that the second (tensor) term in eq. (|30| ) is negligible at large N c . Using the identity 



^ fih ^ 92 =&i-^x^ A y 9 t (33) 

one can rewrite the leading (first) term of eq. (0) as 



+ (^t 75 ^ )2 _ { ^ T A^2 _ ^^^2 (34) 

which resembles closely the Nambu- Jona-Lasinio model. It should be stressed though that in contrast 
to that at hoc model the interaction ( |3"0"| ) i) violates explicitly the Ua(1) symmetry, ii) has a fixed 
interaction strength and Hi) contains an intrinsic ultraviolet cutoff due to the formfactor function 



M(k). This model is known to lead to chiral symmetry breaking, at least at large N c when the use 
of the mean field approximation to the model is theoretically justified. 

The bosonization of these interactions has been performed in (Diakonov and Petrov, 1986a); it paves 
the way to studying analytically various correlation functions in the instanton vacuum. 

A separate issue is the application of these ideas to hadrons made of light and heavy quarks (Cherny- 
shev, Nowak and Zahed, 1994). 



all degrees of freedom in QCD can be divided into two categories: i) those with masses > 1/p and 
ii) those with masses < If one restricts oneself to low-energy strong interactions such that 

momenta are <C 1/p — 600 MeV, one can neglect the former and concentrate on the latter. There 
are just two kind of degrees of freedom whose mass is much less than the inverse average size of 
instantons: the (pseudo) Goldstone pseudoscalar mesons and the quarks themselves which obtain a 
dynamically-generated mass M ~ (l/p)(p 2 /R 2 ) <C 1/p. Thus in the domain of momenta k <C 1/p 
QCD reduces to a remarkably simple though nontrivial theory of massive quarks interacting with 
massless or nearly massless pions. It is given by the partition function (Diakonov and Petrov, 1984b, 



Notice that there is no kinetic energy term for the pions, and that the theory is not a renormalizable 
one. The last circumstance is due to the fact that it is an effective low-energy theory; the ultraviolet 
cutoff is actually 1/p. 

There is a close analogy with solid state physics here. The microscopic theory of solid states is 
QED: it manages to break spontaneously the translational symmetry, so that a Goldstone excitation 
emerges, called the phonon. Electrons obtain a "dynamical mass" m* due to hopping from one atom 
in a lattice to another. The "low energy" limit of solid state physics is described by interactions of 
dressed electrons with Goldstone phonons. These interactions are more or less fixed by symmetry 
considerations apart from a few constants which can be deduced from experiments or calculated 
approximately from the underlying QED. Little is left of the complicated dynamics at the atom 
scale. 

What Petrov and I have attempted, is a similar path: one starts with the fundamental QCD, finds 
that instantons stabilize at a relatively low density and that they break chiral symmetry; what is left 
at low momenta are just the dynamically massive quarks and massless pions. One needs two scales 
to describe strong interactions at low momenta: the ultra-violet cutoff, whose role is played by the 
inverse instanton size, and the dynamical quark mass proportional to the square of the instanton 



10 QCD at still lower energies 



Using the packing fraction of instantons p/i?~l/3asa new algebraic parameter one observes that 



1986a) 




(35) 



density. If one does not believe our variational calculations of these quantities one can take them 
from experiment. 

If one integrates off the quark fields in eq. fl35|) one gets the effective chiral lagrangian: 



7 M, 



S eff [ir} = — iV c Tr In (i ft + iMW 5 
U = exp(i7r A r A ), £P 5 = exp(z7r A r A 7 5 ), L„ = iU ] d^U. (36) 

One can expand eq. (|36p in powers of the derivatives of the pion field and get: 

F 2 r N r 

Sefffr] = ^Jd 4 x TiLl - J d*x [2Tr(^L M ) 2 + TrL^L^ 

+ 240^2 J d x e <xPi5t TrLaLpLjLgLe + ... (37) 
The first term here is the old Weinberg chiral lagrangian with 

2 r d 4 k M 2 (k) . . 

Fl = ±N c —-- ■ t \> - (38) 



(27r) 4 [A; 2 + M 2 (fc)] 2 ' 

the second term are the four- derivative Gasser-Leutwyler terms (with coefficients which turn out 
to agree with those following from the analysis of the data); the last term in eq. (p7|) is the so- 
called Wess-Zumino term. Note that the F n constant diverges logarithmically at large momenta; the 
integral is cut by the momentum- dependent mass at k ~ 1/p, so that one gets the same expression 
as in a different approach described in section 10, see eq. (EH). 



An ideal field of application of the low-momentum partition function (|35|) is the quark-soliton 
model of nucleons (Diakonov and Petrov, 1986b) - actually the model has been derived from this 
partition function. The size of the nucleon ~ (250 MeV)^ 1 is much larger than the size of instantons 
~ (600 MeV) -1 ; hence the low-momentum theory (|35|) seems to be justified. Indeed, the computed 
static characteristics of baryons like formfactors, magnetic moments, etc., are in a good accordance 
with the data (see the talk by K.Goeke). 



11 Can instantons give confinement? 



Our analytical calculations sketched in these lectures, the extensive numerical studies of the instanton 
vacuum by Shuryak and collaborators and the recent direct lattice measurements - all point out that 
instantons play a crucial role in determining the world of light hadrons, including the nucleon. 
Confinement has not much to do with it - contrary to what has been a common wisdom a decade 
ago and in what many people still believe. Nevertheless, confinement is a property of QCD, and 
one needs to understand the confinement mechanism. What can be said today is that confinement 
must be " soft" : it should destroy neither the successes of the perturbative description of high-energy 
processes (no "string effects" there) nor the successes of instantons at low momenta. 



Quite recently we have noticed (Diakonov and Petrov, 1995) that one can obtain the area behaviour 
of large Wilson loops (i.e. confinement) from instantons, provided their size distribution behaves 
as ~ dp/p 3 for large p. It is remarkable that only the large-distance behaviour of the interquark 
potential is sensitive to the tail of the size distribution - all other quantities discussed in this paper 
are determined rather by the average sizes. Moreover, precisely the "one over cube" distribution 
law seems to follow from the statistical mechanics of instantons, when one takes into account the 
quantum gluon fluctuations in the instanton ensemble. 

My special gratitude is to Victor Petrov with whom we worked together for many years on the topics 
discussed in this talk. 
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